Quantum information processing via a lossy bus 
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We describe a method to perform two qubit measurements and logic operations on pairs of qubits 
which each interact with a harmonic oscillator degree of freedom (the bus), but do not directly 
interact with one another. Our scheme uses only weak interactions between the qubit and the bus, 
homodyne measurements, and single qubit operations. In contrast to earlier schemes, the technique 
presented here is extremely robust to photon loss in the bus mode, and can function with high 
fidelity even when the rate of photon loss is comparable to the strength of the qubit-bus coupling. 

PACS numbers: 03.67.Lx, 42.50.Pq, 85.25.Cp, 78.67.Hc 



The practical implementation of quantum information 
processing (QIP) devices is one of the central aims of 
quantum information science. Large scale quantum com- 
puters can dramatically outperform classical computers 
for certain problems while few qubit devices would 
allow other useful QIP applications, including quantum 
repeaters and cryptography 

Of particular interest are bus mediated QIP architec- 
tures, in which coherent quantum operations on multiple 
qubits are implemented indirectly, via a continuous de- 
gree of freedom (the bus) such as a quantized electromag- 
netic field mode which interacts with each of the qubits. 
Such architectures have a number of potential advantages 
for practical implementations. Firstly, in some systems, 
engineering a coherent interaction between qubits and a 
bus mode can be easier than directly coupling pairs of 
qubits. Secondly, architectures involving a bus mode can 
be more readily scaled to large numbers of qubits. When 
a new qubit is added to the system, only it's interaction 
with the bus mode needs to be calibrated and controlled. 
Bus mediated operations can then be implemented on 
any pair of qubits, rather than being restricted to (say) 
nearest neighbour qubits. Finally, in many cases the bus 
can be propagated long distances, which can be useful 
for distributed applications such as quantum repeaters. 

Recently, additional interest in bus-mediated QIP 
has been generated by several experiments which have 
demonstrated the coherent coupling between qubits and 
continuous degrees of freedom. These experiments in- 
clude the strong coupling of a superconducting charge 
qubits to a single electromagnetic field mode the 
coupling of trapped ion qubits to a collective vibrational 
mode [5] and the observation of strong couplingin various 
quantum dot-cavity and atom-cavity systems pi, 0, fill- 

Previously, one of us introduced a scheme for imple- 
menting two-qubit operations via a bus mode fioj . This 
was described in terms of an all-optical implementation, 
using a weak cross-Kerr interaction between the photonic 
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qubits and a bus which is a single mode of the optical 
field, though the scheme could be modified to work for 
other systems with a suitable qubit-bus coupling Hamil- 
tonian. In this scheme, the bus mode is initially pre- 
pared in a coherent state |a) = e~' a ' I 2 ^2 n a n \n) /Vri\, 
where \n) denotes the n'th excited state of the oscillator. 
Subsequently the bus interacts with the qubit degrees 
of freedom via appropriately chosen cross-Kerr interac- 
tions. Finally, by performing a homodyne measurement 
on the bus degree of freedom, the bus is disentangled from 
the qubits, and the entanglement is transferred to the 
qubits. In the original proposal [lfj, this technique was 
used to implement two-qubit entangling measurements 
on the qubits, which arc sufficient for universal quantum 
computation [llj]. A modification of the scheme also al- 
lows the implementation of a CNOT gate [T3|. Other 
schemes for implementing gates using a bus mode have 
also been proposed [H, [Tj| . 

However, these schemes have a substantial drawback 
with regard to practical implementations. At various 
stages in the protocol, the bus mode is in a superposi- 
tion of coherent states, \ae l0q ), with distinct phases {9 q } 
which correspond to different states of the qubits. Super- 
positions of this form are rather fragile in the presence of 
dissipation (e.g. photon loss) of the bus degree of free- 
dom. Hence, losing even a small number of quanta from 
the bus can lead to complete decoherence of the combined 
qubit-bus state, and destroy any entanglement between 
the qubits. A careful analysis shows that the scheme in- 
troduced in , and those derived from it [l2j , can only 
work in the regime of extremely low dissipation defined 
by k/x ^ X t j where k is the loss rate for quanta in 
the bus mode, x 1S the interaction strength between the 
qubits and the bus, and r is the corresponding interac- 
tion time. Since, in most realistic implementations, the 
bus mode will be subject to dissipation, this requirement 
can be a significant restriction. 

In this Letter, we describe a method for bus mediated 
quantum information processing which is robust to dis- 
sipation of the bus mode. The physical resources used 
for this scheme are similar to those used in the original 
scheme namely preparation of a coherent state of 
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the bus mode, conditional phase shifts of the bus and 
homodyne measurements. In addition, the new scheme 
also makes use of unconditional displacements of the bus 
mode; such displacements are straightforward to imple- 
ment in most systems, e.g. by driving the bus mode with 
a classical field resonant with the oscillator frequency. In 
particular, we show how, with these resources, a near- 
deterministic, two qubit projective parity measurement 
(i.e. a measurement of the operator Z\Z<i) can be im- 
plemented. Such a measurement, when augmented with 
single qubit operations, is universal for quantum compu- 
tation, either using the CNOT construction of Ref. [15| 
or by constructing cluster states [H, fl7f . The interac- 
tion between the qubits and the bus can be controlled in 
such a way that the scheme is robust to losses in the bus 
mode even in the regime k > x, where the strength of 
the qubit-bus interaction is weak compared the loss rate. 
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FIG. 1: The photon loss tolerant scheme, (a) The 'circuit' 
that implements the scheme, (b) A phase space illustration of 
a single round of the protocol. In the absence of dissipation, 
we would have 8' — —8, /3i2 = —2aiCOs8 and P34 = 2a?;. 
Dissipation leads to a slight deformation of the bus' trajectory. 

Our scheme for implementing parity measurements is 
illustrated in Fig. [TJ The bus is prepared in the coherent 
state \oti), with cti real. The bus then interacts with the 
qubits via repeated applications of the circuit shown in 
Fig. [BJa). The Uj(9) operation is a conditional phase 
shift operation, which acts on qubit j and the bus mode 
as |0)|a) -> \0)\ae l9 ) and |l)|a) -> \l)\ae' l(> ) (through- 
out this Letter, the states are written with the ket for 
the bus mode at the right). D{(3) = e' 3a+_ ' 3 is the dis- 
placement operator acting on the bus mode. The effect 
of a single round of the protocol can be understood by 
considering the phase space diagram in Fig. [DJb). For 
the 2-qubit input states 1 01) and |10), the bus mode fol- 
lows the indicated loops in phase space. 6, 9', (3 12 and 
/?34 are chosen such that the bus returns its initial state, 
\ati), at the end of each round of the protocol. For the 
input states 1 00) and |11), the consecutive Uj(9) opera- 
tions cancel each other out, and the bus mode follows a 



path along the x— axis. Since ^34 > (3\2, the final state 
of the bus, for the even parity states, is displaced by a 
small distance 8 ps oti9 2 . By repeating this procedure 
N times, the relative displacement of the bus mode for 
the odd and even parity qubit inputs can be made suffi- 
ciently large that the two cases can be distinguished by a 
homodyne measurement of the x— quadrature of the bus 
mode. Finally, if the measurement outcome indicates an 
odd parity state, a fixed rotation about the z— axis is 
applied to one of the qubits to correct for an unwanted 
phase accumulation during the protocol. 

As we show below, by making 9 sufficiently small, this 
procedure can be made extremely tolerant to photon loss. 
Heuristically, this works because, for a given input state 
parity, the bus mode is never in a superposition of co- 
herent states with substantially different complex ampli- 
tudes. Thus, while a moderate amount of photon loss 
reduces the total amplitude of the bus state, it does not 
destroy the coherence of the superposed states within the 
relevant parity subspaces. 

To analyze this scheme in more detail, we consider a 
pair of qubits coupled to a harmonic oscillator (the bus) 
via the interaction picture Hamiltonian (h = 1) 
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Here, Xi(t) denotes the strength of the coupling between 
qubit i and the bus mode, Zi is the Pauli z— operator for 
the ith qubit, a is a lowering operator for the bus, and 
e(t) denotes the strength of the external driving field act- 
ing on the bus, which is resonant with the bare oscillator 
frequency [2lJ. We assume the phase and amplitude of 
this driving term can be accurately controlled. The first 
two terms in Eq. (fTJ) generate the conditional phase shifts 
of the oscillator, while the third term generates the un- 
conditional displacements of the oscillator. This Hamil- 
tonian can be realized in a variety of settings, including 
the dispersive regime of cavity quantum electrodynamics 
(CQED) [22J (which can be realized in a variety of atomic 
and solid-state systems, and can also be simulated in an 
ion trap system [al), and in an all-optical setting, via the 
cross-Kerr effect (lQj ] . 

A single round of the scheme can be implemented 
by a piecewise constant time dependence for Xi(t) an d 
e\t). We take xi(*) = -X2CO = {Xo, 0, ~Xo, 0}, and 
e(t) = {0, — £0, 0, £0} for the intervals {0 < t < t\, t\ < 
t < i 2 , h < t < is, t 3 < t < i 4 }. This has the effect 
of implementing the consecutive pairs of Uj(±9) opera- 
tions [shown in Fig. HJa)] simultaneously, which is con- 
ceptually equivalent to performing them consecutively, 
since the Hamiltonian terms that generate these oper- 
ations commute with one another. We also define the 
time intervals Tj = t j — corresponding to the length 
of each step of the evolution. The results obtained using 
this choice for Xi{t) and e{t) will be qualitatively valid 
for other choices, too. For example, in a particular im- 
plementation, it may be easier to have a fixed, constant 
qubit-bus coupling, and to periodically apply very strong, 
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short driving pulses (with |e| 3> \Xi\) to the bus mode. 

We describe dissipation using a standard approach for 
the CQED setting: the damping from the oscillator may 
be modeled as leakage from the corresponding cavity mir- 
ror at a rate 2k. Within the Born, Markov, and rotating 
wave approximations, such losses can be described by the 
non-unitary term in a Lindblad master equation [181 ]: 



-2kt 3 \ 



p(t) = -i [Hi(t), pit)} + 2KD[a]p(t), 



(2) 



where p(t) is the density matrix of the combined system 
of both qubits and the oscillator mode, and £>L4]/c> = 
ApA^ - (pA^A + A*Ap)/2 Below, we solve Eq. © 
for a single round of the protocol, and then use the results 
to determine both the total number of rounds required to 
distinguish even parity states from odd parity states, and 
the total dephasing error at the end of these N rounds. 

The solution of Eq.© can be simplified by noting that, 
since the aim of the scheme is to perform a parity mea- 
surement, only dephasing errors within the correspond- 
ing parity subspaces are of concern. Thus we are only 
interested in solutions of Eq. © for initial states of the 
form O|00) +y\ll))\a l ) and (a;|10) + y\01))\a l ), where 
x and y are arbitrary amplitudes. We take the worst 
case, x = y = l/\/2, and so the resulting errors should 
be viewed as an upper bound. The solution can be fur- 
ther simplified by projecting p(t) onto the computational 
basis of the qubits. This results in a set of uncoupled dif- 
ferential equations for the operators Pi\j(t) — (i\p(t)\j), 
where i,j £ {00, 01, 10, 11}, which can be solved indepen- 
dently. These equations are analogous to those describ- 
ing the single-qubit, driven Jaynes-Cummings model in 
the dispersive limit, and thus can be solved exactly us- 
ing techniques similar to those used in Refs. 0, [20j | . 
We omit the details of this calculation, but give the full 
analytic results here. 

For an input state |\l/ e+ )|ai), where |* e ±) = (|00) ± 
\ll))/y/2, the state after a single round of the protocol 
(i.e. at time (4), is given by |\I/ e +)|a e ), with 
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where (5\ 2 = — (co/k)[1 - e 
and the time intervals t» = t< — U-i correspond to the 
length of each step of the evolution. 

For an input state |\& 0+ )|aj), where \^ ±) = ( 1 01) ± 
|10))/a/2, the state after a single round is 
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1)], r 12 = a 2 e" 2KTl (l - e~ 2KT2 )(l - e 2 ^), $i 2 
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corresponds to a dephasing er- 



e -K(ri+T 2 ) e ixTl Q , i 

— Toi — Ti 2 — $i 2 — r 2 3 
ror (see below), while the imaginary part leads to a de- 
terministic phase rotation which is ultimately corrected 
by applying the single qubit unitary U z [see Fig. Q] (a)], 
after the homodyne measurement stage of the protocol. 

The timcscalcs t 2 , r 3 and T4 are chosen as follows. r 2 
and T3 can be fixed requiring that, for symmetry, at time 
£3, the bus mode associated with the input state |\? ±) 
is in the state a 3 = —on. This leads to the approximate 
values t 2 w (2 — k 2 t 2 — x 2 T 2 )ai/eo and t 3 w n(l — 2kt\). 
T4 is chosen such that a = on, i.e. that the bus mode 
associated with the input state | v E'o±) is returned to it's 
initial state at the end of each round of the protocol. 

Using these values of tj, it is possible to obtain 
approximate expressions for the total dephasing error 
at the very end of the protocol. After N rounds, 
the state corresponding to an initial state \^o+) is 

pi N) = [(i-p)|f +)(f 0+ |+p|* -)(* -l] ® Wi){ca\, 

where p = (1 - e" ArRc ( r ° 1+ri2+r23 ))/2 is the probability 
of a dephasing error, and |* ±) corresponds to |^ ±)j U P 
to the single qubit correction U z discussed above. The fi- 
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Writing the expression for the dephasing error in 
terms of t\ gives the approximate expression p s» 
4A r a 2 «x 2 ' r i/3 + 4iVafKX 2 T 2 /e - The first term in this 
expression originates from decoherence on the 'rotation' 
parts of the evolution (i.e. to — > t\ and t 2 — > t 3 ), whereas 
the second term originates from the first 'displacement' 
part [t\ — * t 2 ). Note that no decoherence results from 
the second displacement stage (i 3 — > t±), as the bus mode 
is not in a superposition state during these parts of the 
evolution. 

The total number of required rounds of the protocol, 
N, is fixed by the requirement that |q4 ) and |aj) are 
distinguishable by a X— homodyne measurement at the 
end of the protocol. For coherent states, this implies that 
ote^ — oti — A, where A > 1 is a constant that depends on 
the required error in distinguishing the states |^ ±) from 
the states |\P e ±)< This requirement gives N = A/a^x 2 ^- 
Thus the total error is well approximated by 
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(4) 



The first term in this expression can be made arbitrarily 
small by decreasing the value of kt\ . The second term in 
this expression is proportional the ratio ft/eo- Thus errors 
can be minimized provided the strength of the classical 
driving field can be made large compared to the pho- 
ton loss rate. An important point to note is that both 
terms are independent of the ratio k/x> an d therefore 
the scheme is robust to losses in the bus mode even in 
the regime k > x, where the strength of the qubit-bus 
interaction is weak compared to the loss rate. The ap- 
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proximate error, together with the exact values for the 
same values of the parameters r i; are plotted in Fig. 
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FIG. 2: (a) Exact (curves) and approximate (points) values 
of the final dephasing error against ri. Other parameters: 
A = 1.5, on = 2, k = 0.05, x = l, e = 10 2 ,10 3 ,10 4 , (top 
to bottom), (b) Overhead, N, required to reach a given error 
rate, p, for A = 1.5, on = 2, k = 0.05, x = 1, e = 10 5 . 

The price to be paid for this robustness is an increase 
in the overhead cost of the scheme: the first term in Eq. 

can be reduced by reducing n, but as t\ is reduced, 
the required number of rounds of the protocol required 
increases. If p is dominated by the first term, then the 
required number of rounds is N (16 A 3 at / '9) (k / x) 2 /p 2 - 
Thus the overhead cost grows quadratically with the re- 
quired reduction in the errors (see Fig. 

In view of this overhead, it is worth briefly consider- 
ing the effect of errors in the scheme. Stochastic noise in 
the control pulses used to implement the protocol tend 
to broaden final state of the bus mode, making the odd 
and even parity outcomes harder to distinguish by the 
X— homodyne measurement. Treating the noise in each 
pulse as independent Gaussian processes, the final uncer- 
tainty in the real part of is found to be a 2 (aQ N ^) = 



Na 2 (a ) « N{a 2 ((3 12 ) + a 2 {(3 M ) + (3 2 l2 sin 2 0a 2 {6)\ » 
2Na 2 (/3 12 ) +4a Aa 2 (9), where a 2 (. . .) denotes the vari- 
ance in each parameter. Thus the noise will not signif- 
icantly affect the outcome of the X-homodyne measure- 
ment provided o-((3i 2 ) ^ and a{9) < (4ao^4)~ 5 . If 
these conditions cannot be met then the scheme can still 
be made to work by increasing the number of rounds, 
since the peak separation grows linearly in N, while the 
peak width grows only as . Another potentially im- 
portant source of noise is single qubit errors. These can 
potentially build up during the implementation of the 
two-qubit operation, leading to an effective error rate 
that scales with N. Thus our scheme is particularly 
suited to systems in which the single qubit decoherence 
rate is small compared to the bus damping rate. Fortu- 
nately, there are many systems (such as trapped ions [5j]) 
where the single qubit error rate is much longer than all 
other timescales. 

In conclusion, we have proposed a method for imple- 
menting two-qubit operations via a quantum bus, which 
is extremely robust to dissipation of the bus degree of 
freedom. Our method implements a parity measurement 
of the two-qubit system, which, when augmented with 
single qubit operations, is sufficient for universal compu- 
tation. The residual errors in our scheme can be made 
negligible even when the strength of the qubit-bus inter- 
action is weak compared to the loss rate. We anticipate 
that this robustness will permit bus-mediated QIP in a 
variety of physical setups, including superconducting sys- 
tems, trapped ions or atoms, and all-optical systems. 
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